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We calculate parity-dependent level densities for the even-even isotopes 58 ' 62 > 66 p e anc j 58 Ni and the 
odd-A nuclei 59 Ni and 65 Fe using the Shell Model Monte Carlo method. We perform these calcu- 
lations in the complete fp-gds shell-model space using a pairing+quadrupole residual interaction. 
We find that, due to pairing of identical nucleons, the low-energy spectrum is dominated by positive 
parity states. Although these pairs break at around the same excitation energy in all nuclei, the 
energy dependence of the ratio of negative-to-positive parity level densities depends strongly on the 
particular nucleus of interest. We find equilibration of both parities at noticeably lower excitation 
energies for the odd-A nuclei 59 Ni and 65 Fe than for the neighboring even-even nuclei 58 Ni and 66 Fe. 



I. INTRODUCTION AND MOTIVATION 

Nuclear level densities play an important role in theo- 
retical estimates of nuclear reaction rates needed in var- 
ious applications including astrophysical nucleosynthesis 
processes like the s-, r-, and rp-process |l| They also con- 
tribute to one of the largest uncertainties in the cross 
section determinations used for large-scale nucleosynthe- 
sis networks [T], @, y, |j| • Typically, nuclear level densities 
are described in these astrophysical studies by using the 
backshifted Fermi gas model of Gilbert and Cameron Q . 
This model extends the non-interacting Fermi gas model 
of Bcthe [6| by considering pairing among like nucleons 
via a backshift of the excitation energy E. This backshift 
accounts for the energy required to break nucleoli pairs. 
Furthermore, the astrophysical applications assume an 
empirical ansatz for the angular momentum distribution 
of the levels and consider equilibration of both parities 
at the energies of interest [J, |j| . 

For astrophysical applications, one is often interested 
in the nuclear level density at rather low excitation en- 
ergies. For example, the typical neutron energies in 
r-process nucleosynthesis reach only up to a few MeV 
energies, nuclear structure and pairing 
effects strongly influence the level density and an equili- 
bration of both parities is quite unlikely. In particular, in 
even-even nuclei, pairing among identical nucleons gener- 
ates only positive parity states even when single-particle 
states of opposite parity are present. Therefore, negative 
parity states should be suppressed in even-even nuclei 
at low energies due to both pairing and the underlying 
single-particle structure of the mean field which groups 
states of the same parity, at least for nuclei with mass 
numbers A smaller than about 100. The parity equili- 
bration is governed by the energy scales associated with 
pair-breaking and the shell gap between opposite-parity 
states near the Fermi surface. While the latter is typi- 
cally of order 5-6 MeV for intermediate mass nuclei, the 
former strongly depends on the nuclear structure. 



We investigate in this paper the competition of pairing 
and single-particle structure on the parity dependence of 
the level density. We perform Shell Model Monte Carlo 
calculations of the even-even nuclei 58 > 62 . 66 F e anc j 
and the odd-A nuclei 59 Ni and 65 Fe in the complete fp- 
gds model space, where the single-particle states in the 
fp shell have negative parity, while those of the gds shell 
have positive parity. As we shall see, residual interactions 
among the protons and neutrons also influence the level 
density. 



II. FORMALISM 

The Shell Model Monte Carlo (SMMC) approach al- 
lows the determination of nuclear properties at finite tem- 
perature in unprecedentedly large model spaces, consid- 
ering the important correlations among the valence nucle- 
ons @, H, The SMMC method describes the nucleus 
by a canonical ensemble at temperature T = 0~ 1 and 
employs a Hubbard-Stratonovich linearization [lCj of the 
imaginary-time many-body propagator, e~@ H , to express 
observables as path integrals of one-body propagators in 
fluctuating auxiliary fields 0, [H|- Since Monte Carlo 
techniques avoid an explicit enumeration of the many- 
body states, they can be used in model spaces far larger 
than those accessible to conventional methods. The no- 
torious sign problem that plagues Monte Carlo studies of 
Fermionic systems can be avoided within the SMMC by 
adopting a pairing+quadrupole force as residual interac- 
tion [8j. 

There are already several successful studies of level 
densities using the SMMC approach. Ormand [T^, as 
well as Nakada and Alhassid jl3j |. calculated the level 
density for selected even-even nuclei, while Ref. [3] ex- 
tended this work to odd-A and odd-odd nuclei. Using 
projection techniques, Alhassid and collaborators suc- 
ceeded in studying the parity dependence fl3j and the 
angular-momentum dependence [161 ] of the level density 
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for intermediate-mass nuclei. Recently, Ref. [ljj ex- 
plored the influence of pairing correlations and their en- 
ergy dependence on nuclear level densities. All these ap- 
proaches are based on the ability of the SMMC to calcu- 
late expectation values of an observable at temperature 
T as the thermal average of a canonical ensemble (with 
fixed numbers of protons and neutrons). Thus for the 
energy excitation function, one has 



E(p) = 



Tr [He-? H ] _ Tr [He^ H ] 



Tr[e- 



Z((3) 



(1) 



where Z is the partition function, (3 the inverse temper- 
ature, and H the nuclear Hamiltonian. Using the parity 
projection operators P± = (1 ± P)/2, where P is the 
parity operator, one is able to calculate the ratio of the 
parity-projected partition function to the total partition 
function fl3j 



Z±{fi) Tr[P±e-P B ] 



Z(f3) 



(2) 



from which one can then extract the parity-projected en- 
ergy excitation functions 



E±(p) 



E{f3) 



d\nY±(f3) 
dp 

Tr [HP±e~l 3H ] 

zIW) 

JdE'e-P E 'Ep±(E') 



(3) 



where we have introduced the parity-projected partition 
functions Z± and level densities p±. 

To obtain p± from Eq. [3] requires an inverse Laplace 
transform which we treat within the saddle-point approx- 
imation: 



P±(E) = 



a /3E ± +lnZ ± (l3) 



-2tt 



dE ± {!3) 
d/3 



(4) 



where (3 = (3{E±) is obtained by inverting E±({3). An 
analog relation holds between the total level density p(E) 
and the energy excitation function E((5). 



III. RESULTS 

Our SMMC calculations were performed in the com- 
plete fp-gds model space with 50 valence orbitals for 
both protons and neutrons. The single-particle en- 
ergies and the residual pairing+quadrupole interaction 
were adopted from previous SMMC studies that success- 
fully explored the competition of isovector pairing ver- 
sus quadrupole deformation in the A ~ 80 mass region 
[HI], and as a function of temperature in selected nu- 
clei. These previous calculations clearly identified the 
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FIG. 1: Parity-projected level densities for 58 > 62 - 66 Fe and 58 Ni. 



breaking of pairs of identical nucleons in even-even nu- 
clei around T w 0.7 MeV 0. 

In the following, we discuss various results of these cal- 
culations. Figure Q] shows the total and parity-projected 
SMMC level densities for 58 . 62 > 66 Fe and 58 Ni. We observe 
that, at modest excitation energies, the total level density 
increases with the mass number. This is expected from a 
Fermi gas approximation to the level density, which scales 
like p ~ exp{2\/ aE} with the level density parameter a 
proportional to the mass number A [l[ . Interestingly, we 
observe an equilibration of both parities at E as 11-12 
in 58 Fe, while equilibration occurs already around £«6 
MeV in 62,66 Fe. This observation simply reflects the fact 
that it is energetically more costly to promote neutrons to 
the gds shell for 58 Fe, which requires neutrons to bridge 
the gap from the p 3 / 2 to the g 9 / 2 orbital, than for 62,66 Fe. 
Nevertheless, it is not simply the occupation number of 
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the (?g/2 orbital which matters. This quantity is shown in 
Fig. [1 For 58 Fe, the SMMC predicts rather little aver- 
age occupation of the g 9 / 2 orbital, which increases from 
0.2 neutrons in the ground state to 0.3 at T = 1 MeV 
(corresponding to an excitation energy of about E ~ 6 
MeV) to 1.0 at T = 2 MeV (E ~ 23 MeV). On the other 
hand, one finds an occupation number of order 1.0 in the 
62 Fe ground state, which grows to 2.1 at T — 2 MeV. For 
66 Fe, the residual interaction promotes neutrons to the 
gg/2 orbital, yielding a g 9 / 2 occupation number of about 
3.0. The occupation grows only mildly with tempera- 
ture and reaches about 3.6 at T = 2 MeV. We note that 
the residual interaction mainly scatters (neutron) pairs 
from the fp shell to the gds shell. Although these corre- 
lations are important for the nuclear structure, they do 
not change the parity. This is clearly seen in a compar- 
ison between 62 Fe and 66 Fe, which shows quite a similar 
ratio of negative-to-positive level densities (Fig. []}, while 
clearly these nuclei have distinct occupation numbers of 
the <?g/2 orbital (Fig. [2]). While a promotion of nucleons 
(neutrons) is required to make negative-parity states in 
our model space, the strong pairing correlations also have 
to be overcome. 

As in [ijj [2(|, we define the pairing strength as the 
sum over all matrix elements of the pair matrix 



{A\j a ,j b )A{j c ,j d )), 



with the J = pair operator 
1 



a i x a i 



JM=00 



(5) 



(6) 



where aj a creates a nucleon in the orbital a with angular 
momentum j a [22| . Since only genuine pair correlations 
are of interest, we subtract the pure 'mean-field' values, 
i.e. those obtained without residual interaction. 

As is shown in Fig. [3J the pairing strengths decrease 
strongly with temperature, and are reduced to about half 
of the ground-state values around T — 0.8 MeV (corre- 
sponding to E 4 MeV). This breaking of pairs is ac- 
companied by a peak in the specific heat (e.g. [20j]). We 
note that the proton pairing strengths are quite similar 
for all three iron isotopes (e.g. [i~4j). 

For 62 Fe and 66 Fe the ratio of negative-to-positive level 
densities is quite similar and equilibration is reached at 
excitation energies around 6 MeV; i.e., once the pairs 
with positive parities are sufficiently broken. This ap- 
parently requires more energy than the gap between the 
last occupied neutron orbital and the <7g/ 2 orbital with 
its large degeneracy. Hence, once pairs are broken, many 
negative-parity states can be formed in 62 Fe and 66 Fe. 
The situation is obviously different for 58 Fe where the 
gaps from the last occupied proton (/V/2) and neutron 
{Pz/2) orbitals to the g 9 / 2 orbital are larger than 6 MeV. 
Thus, even at the energy at which pairs are broken, it is 
more likely to form positive-parity states with unpaired 
nucleons in the fp shell than negative-parity cross-shell 
states. 
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FIG. 2: Average SMMC occupation numbers of the gg/2 neu- 
tron orbitals for 58 - 62 > 66 F e as function of temperature. 
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FIG. 3: Pairing strength for neutrons in 58 > 62 ' 66 F e as function 
of temperature. As explained in the text, the pairing strength 
is defined as the sum of the matrix elements of the pair matrix 
[5] with the 'mean-field' value subtracted. 



Ref. [IH describes a model to estimate the ratio of 
negative-to-positive level densities on the basis of the in- 
dependent particle model and a BCS treatment of pair- 
ing. It appears that this model slightly underestimates 
this ratio. This is due to the fact that the residual inter- 
action mixes orbitals from adjacent shells with different 
parities at lower energies than obtained in the indepen- 
dent particle model. 

Although they have the same neutron number, the iso- 
tones 58 Fe and 58 Ni differ in their single-particle struc- 
ture, and this shows in their level densities. The fact 
that 58 Fe has four neutrons in the P3/2 orbital, compared 
to two in the case of 58 Ni, allows for a larger amount of 
negative-parity states at low energies in 58 Fe, since these 
states are generated by configurations with one and three 
neutrons in the gds-orbitals for 58 Fe but with only one 
neutron for 58 Ni. As a consequence, the parity unbalance 
in the level density persists to higher excitation energies 
in 58 Ni (about E x = 16 MeV, see Fig. P than in 58 Fe. If 
the total level densities of the two nuclei are compared, 
the one for 58 Fe is noticeably larger at a given excitation 
energy than the one for 58 Ni. At modest excitation ener- 
gies, this difference can be approximately accounted for 
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by a constant energy shift of about 3 MeV, simulating 
the shell gap in 58 Ni. 

While the structure of even-even nuclei is strongly 
dominated by pairing at low excitation energies, thus 
leading to the dominance of positive-parity over negative- 
parity states, the situation is different for odd-^4 nuclei. 
In this case, the unpaired nucleon is not hindered by 
pairing and can, depending on the single-particle struc- 
ture, occupy negative-parity states (fp shell) or positive- 
parity states (gds shell). To investigate the differences 
between even-even and odd-A nuclei, we also calculated 
the SMMC level densities for 59 Ni and 65 Fe. For 59 Ni, 
the unpaired neutron occupies ap 3 / 2 orbital. Despite the 
fact that this neutron is not hindered by pairing, it will 
nevertheless, at low excitations energy, mainly occupy or- 
bitals in the fp shell; consequently, negative-parity states 
dominate the 59 Ni level density at low excitation energies. 
This is indeed born out of our SMMC calculation (Fig. 

which yields a balance between negative- and positive- 
parity level densities at excitation energies around 13 
MeV. Comparing with 58 Ni (Fig. [T|), where the equili- 
bration energy is about 16 MeV, we see that neutron 
pairing, e.g., the residual interaction, indeed has a sig- 
nifant effect in the equilibration of positive and negative 
parities. 

The neutron single-particle structure of 65 Fe corre- 
sponds to a single hole in the fp shell. Thus, low-energy 
excitations can be achieved by either changing the hole 
structure in the fp shell or by promoting the neutron 
to the gg/2 orbital, which, however, has the opposite 
(positive) parity to the ground state. As a consequence, 
one expects equal distribution of negative- and positive- 
parity level densities already at low excitation energies. 
This is indeed observed in our SMMC calculations (Fig. 
@]), where we find the same level densities for both parities 
down to E x w 1.5 MeV, corresponding to the lowest tem- 
perature for which we have been able to perform SMMC 
calculations in the odd-A systems. (An odd-A sign prob- 
lem occurs even in the presence of good sign interactions 

Bo 



IV. CONCLUSIONS 



strongly on the single-particle structure of the nuclei. If 
the Fermi energies of protons and neutrons are relatively 
well separated from orbitals with the opposite parity, an 
equal distribution between parities in the level densities is 
achieved at higher excitation energies than for nuclei with 
Fermi energies close to the oscillator shell closure. As ex- 
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FIG. 4: Parity-projected level densities for 65 Fe and 59 Ni. 



amples for this difference, we presented parity-projected 
level densities for 58 Fe, 58 Ni, and 66 Fe. The same general 
trend can be observed for odd-A nuclei ( 59 Ni and 65 Fe). 
However, as these nuclei have an unpaired neutron, the 
balance between negative- and positive-parity level den- 
sities is achieved at somewhat lower excitation energies 
than in the neighboring even-even nuclei. 

The SMMC approach has again been proven to be a 
powerful tool to microscopically study nuclear level den- 
sities. In the future, it will be used to explore the angular- 
momentum dependence of the level density. First steps 
towards this goal have been presented in [l6|, [5l[ . 



In summary, we calculated parity-projected level den- 
sities for several even-even and odd- A /p-shell nuclei us- 
ing the Shell Model Monte Carlo approach. For even- 
even nuclei, we confirm that the low-energy spectrum 
is dominated by pairing among nucleons, resulting in 
a dominance of the positive-parity level density. Al- 
though the pairs break at excitation energies of a few 
MeV, the contribution of negative-parity states depends 
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